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Abstract — Carrier aggregation (CA) and heterogeneous base 
stations (HetNets) and are two distinct features of next-generation 
cellular networks. CA increases transmission bandwidth by 
aggregating multiple component carriers on the physical layer, 
while small cells in HetNets are vital for data off-loading and can 
significantly improve area spectral efficiency compared to using 
just macrocells. In this paper, we propose a tractable multi- 
band multi-tier model for CA in HetNets, where base station 
positions in each tier follow an independent Poisson point process 
(PPP). Our model incorporates an appropriate notion of load. 
For a typical HetNet consisting of two tiers, macro and small 
cells, and two bands, 800MHz and 2.5GHz band, we find that 
small biasing factor (less than lOdB) yields little gain and the 
maximum biasing gain is about 25% to 40%, attained at large 
biasing factors. Our model also characterizes the impact of band 
deployment. For the typical HetNet described, our study suggests 
that the 800MHz band be deployed both in macro and small cells 
to ensure coverage and balance network load. Our study also 
indicates that deploying the 800MHz band and 2.5GHz band in 
both the macro and small cells is optimal. 
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Fig. 1: Different CA spectrum scenarios: Component carriers 
can be either intra-band contiguous, intra-band 
non-contiguous or inter-band non-contiguous. 
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I. Introduction 

Carrier aggregation (CA) is considered as a key enabler 
for LTE- Advanced which can meet or even exceed the 
IMT-Advanced requirement for large transmission bandwidth 
(40MHz-100MHz) and high peak data rate (500 Mbps in the 
uplink and 1 Gbps in the downlink). In particular, CA enables 
the concurrent utilization of multiple component carriers on 
the physical layer to expand the effective bandwidth ||2)-||4|]. 
The aggregated bandwidth can be as large as 100 MHz, for 
example, by aggregating 5 20MHz component carriers. CA is 
significantly different from existing multi-carrier transmission 
techniques, e.g., orthogonal frequency division multiplexing 
(OFDM). In particular, each component carrier in CA acts 
as a physical channel and has individual transmission modes 
(modulation and coding scheme). The bandwidth of these 
component carriers can vary widely (ranging from 1.4 MHz to 
20 MHz for LTE carriers). Further, CA can be used in different 
spectrum scenarios including intra-band contiguous CA, intra- 
band non-contiguous CA and inter-band non-contiguous CA 
0, which are illustrated in Fig.[T] The propagation character- 
istics of different component carriers may vary significantly, 
e.g., a component carrier in the 800MHz band has very 
different propagation characteristic from a component carrier 
in the 2.5GHz band. Due to the significance and unique char- 
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acteristics of CA, how should cellular networks be designed 
in order to best exploit CA? 

Cellular networks are undergoing a major evolution as 
current cellular networks cannot keep pace with user demand 
through simply deploying more macro base stations (BS) 
ll6l- lfT0l . As a result, attention is being shifted to deploying 
small, inexpensive, low-power nodes overlaid in the current 
macrocells. These low power nodes may include pico and 
femto BSs, as well as distributed antennas. Cellular networks 
with them take on a very heterogeneous characteristic, and are 
often referred to as HetNets. Not surprisingly, random spatial 
point process, particularly homogeneous Poisson point process 
(PPP) for its tractability, has been used to model the locations 
of the various heterogeneous BSs. Using this approach, many 
works (see, e.g., lfTT1 - lfT4ll ) successfully advance the under- 
standing on HetNets as well as the design. However, CA has 
not been considered in ITT1l - lfT4l yet. 

In this paper we propose a tractable multi-band multi-tier 
model for CA in HetNets. The BS locations are modeled by 
PPPs while different path loss exponents are used to capture 
the possible large differences in the propagation characteristics 
of the available bands. Within this model, the user equipment 
(UE) performs cell selection based on the maximum received 
signal power. Two types of CA are considered: multi-flow 
CA and single flow CA, which will be defined in Section llll 
Essentially, multi-flow CA allows UEs to connect to multiple 
BSs simultaneously while using single flow CA UEs can only 
connect to one BS at one time. The proposed model is flexible 
enough to capture the main unique features of CA but yet 
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tractable enough for analysis. 

One major concern about deploying small cells is that they 
have very small coverage areas due to their low transmit 
power. As a result, small cells are often lightly loaded and do 
not accomplish much, while the macrocells are still heavily 
loaded. To alleviate this issue, biasing allows the small cells 
to expand their areas |6), Q, lfl5l . enabling more UEs to 
connect to small cells. It is expected that biasing can improve 
load balancing and thus leads to higher throughput. However, a 
theoretical study on the impact of biasing is challenging. While 
[ 1 4 1 modeled biasing in a single-band HetNet, it assumes 
a fully-loaded HetNet, and so correspondingly biasing does 
not help. In this paper we introduce an appropriate notion 
of load. Thanks to the proposed load-aware model, we are 
able to perform an extensive study on the impact of biasing 
on HetNet's performance and give useful design insights. For 
example, in a two-tier network consisting of macro and small 
cells we find that biasing does make small cells accomplish 
more when the network load is unbalanced, while biasing 
hurts the network performance in a fully loaded network, 
implying that using a fully loaded model to study biasing is 
inappropriate. Moreover, we find that small biasing factor (less 
than lOdB) yields little gain, while the maximum biasing gain 
is about 25% to 40%, attained at large biasing factors (between 
15dB and 30dB). 

Another novel contribution of this paper is to analytically 
study HetNet's performance under different band deployment 
configurations. With the incorporation of band deployment in 
our model, we perform a comprehensive study in a HetNet 
consisting of two tiers (macro and small cells) and two 
bands (800MHz and 2.5GHz band). Our study suggests that 
the 800MHz band should be deployed in both macrocells 
and small cells. Deploying the 800MHz band in macrocells 
can ensure coverage, while deploying the 800MHz band in 
small cells better balances network load and thus makes the 
small cells accomplish more (even at the cost of increased 
interference for cell edge UEs of small cells). Besides, we 
find that splitting the spectrum between different tiers is 
strictly suboptimal. The reason is that there is a capacity 
loss in macrocells by spectrum splitting and this loss can 
only be partially overcome by adding a large amount of 
small cells. Thus, it has been recognized that HetNets should 
allow full spectrum reuse augmented with appropriate intercell 
interference management iflOl . Indeed, our study shows that 
the deployment that both the two tiers use both the two bands 
is optimal in many scenarios. 

The rest of this paper is organized as follows. Section [TT] 
describes the specific system model and defines the two CA 
schemes. We then present theoretical analysis on the multi- 
flow CA and single flow CA in SectionlHll and ll VI respectively. 
Section[V]presents some numerical results, and is followed by 
our conclusions in Section [VI] 

II. System Model 

We consider a general HetNet with K. = {1, K} denoting 
the set of K tiers which may include macrocells, picocells, 
femtocells, and possibly other elements. BSs of different types 



differ in terms of deployment density, spectrum resource, 
transmit power and supported modulation and coding scheme. 
In this section, we describe the key aspects of this model, 
which includes the BS and UE distribution, radio channels, UE 
association schemes and load modeling, and the performance 
metric. 

A. Distributions of BSs and UEs 

Though regular grid based models are popular for cellular 
networks, BSs cannot be perfectly regularly deployed in reality 
and thus some randomness must be considered. In this paper 
we adopt random spatial point process to model the BS 
locations. More specifically, the BS locations are modeled 
as K independent homogeneous PPPs fl6l . The set of PPP 
distributed BSs in tier k is denoted by $fc with density This 
Poisson assumption on BS distribution originally proposed in 
ifTTl has been recently shown to be about as accurate in terms 
of SINR distribution as the hexagonal grid for a representative 
urban cellular network lfl8l and has been extended to HetNets 
in flj] and lfT3l 

The UEs are also assumed to be randomly distributed 
according to an independent PPP of density \( u \ Or equiva- 
lently, given the number of UEs in a certain area, these UEs are 
uniformly and independently located in this area. This uniform 
distribution assumption on the UE positions is popular and 
mainstream and is often utilized in system level simulations 
by industry. 

B. Channel model 

We assume that there are a set of M available bands (e.g., 
800MHz, 2.5GHz, ...) denoted as M = {1,2,...,M}. The 
bandwidth of each band i is denoted by Bi and the path loss 
exponent of band i is denoted by CKj. Here we use different 
path loss exponents for different bands to capture the possible 
large differences in propagation characteristics associated with 
each band's carrier frequency. 

We further assume that each band is small enough to have 
relatively constant path loss across it. In particular, if band i 
is used by tier k, the received power is modeled as 

Pi, k = p t , k Had- a ' (i) 

where d is the length of the radio link, C{ is a constant that 
gives the path loss in band i when d = 1, H is a random 
variable capturing the fading effects of the radio link, and 
Pi^k is the i-th band transmit power used by BSs in tier 
k, assumed to be fixed. Note that Ci depends strongly on 
carrier frequency, e.g. Ci = (^i/An) 2 where /i^ denotes the 
wavelength. For simplicity, we ignore shadowing and consider 
Rayleigh fading only, i.e., H ~ Exp(l). In fact, as shown 
in lfl9l . the randomness of the BS locations actually helps 
to emulate shadowing. In particular, with large enough shad- 
owing variance (e.g., lOdB) the grid networks "converge" to 
Poisson networks. Thus, incorporating shadowing in Poisson 
networks only enlarges the performance difference between 
grid and Poisson networks. This implies that, to have a better 
performance match between grid and Poisson networks, it is 
not necessary to model shadowing in Poisson networks. 



3 



C. User association schemes 

We consider two types of CA. The first type is called multi- 
flow CA, where UEs can be potentially associated with all 
the available tiers simultaneously (but in different bands) and 
can aggregate data using all the available bands. Assuming 
that the typical UE is located at the origin, in band i the UE 
is associated with tier k that provides the maximum biased 
received power, i.e., k = a,rgma,x(ZiPi,i\\Yi \\~ ai : I e fC), 
where Zi denotes the biasing factor of tier I and Yi denotes 
the location of the nearest BS in tier I. Biasing factors Z = 
{Zk : k e JC] are used in HetNets for the purpose of load 
balancing. In particular, adopting larger Zk expands the cell 
range of BSs in tier fc and thus more UEs can connect to tier k. 
Range expansion is important for HetNets. Small cells do not 
accomplish much without it, since they otherwise have very 
small coverage areas. 

The second type is called single flow CA, where UEs can 
be associated with only one of the available tiers at one time, 
i.e., one BS at some tier, though they still can aggregate data 
using all the available bands used by that tier. For this type 
of CA, the typical UE scans over all the tiers and bands and 
connects to the tier that provides the strongest biased received 
power in some band, say i*. Then the UE performs CA with 
respect to this tier. Formally, the UE connects to tier k such 
that (k,i*) = argmax(Z;.Pu||Yi ir Ql : 6/Cx M). 

In summary, multi-flow allows the UE to perform CA across 
all the available bands in (possibly) different tiersQ while 
single flow allows the UE to perform CA only across the 
available bands by the sole selected tier. Fig. |2] demonstrates 
the model in question and shows the difference between multi- 
flow CA and single flow CA. Note that both schemes are 
based on long term measurements which average out the 
instantaneous effect of fading. This assumption is reasonable 
since BS selection in real cellular networks does average out 
short term fading, which in any case is often exploited or 
suppressed with various forms of diversity. For example, in 
most cases of interest, the bandwidth Bi is much larger than 
the coherence bandwidth of the channel, so the per band 
received power tends to the mean. 

D. Load modeling 

We assume that each UE connecting to tier k in band 
i requires a basic share b of the bandwidth resource B^. 
For example, b may represent the basic frequency resource 
allocation unit in LTE networks. This requirement can be 
satisfied if tier k in band i is under-loaded, e.g., Li : k < Bj/b 
where Li t u denotes the mean number of UEs attempting to 
connect to one tier-fc BS in band i. If Li.k > Bi/b, this 
implies that too many UEs attempt to connect to tier k and 
correspondingly tier k becomes fully-loaded. In this case, 
some UEs within the coverage of tier k have to be blocked. 
Or equivalently, assuming all the UEs are admitted by tier k, 
they can only obtain a fraction of the time domain resource. 

'Note that multi-flow UEs can only connect to one tier in each band i. For 
example, if both the macro and small cells deploy band i, then the UE can 
only communicate with at most one of the nearest macro or small BS in band 
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Fig. 2: System model: Multi-flow versus single flow 



In the following, we stick to the former interpretation and will 
derive the admission probability while keeping in mind that the 
admission probability can also be interpreted as the fraction of 
the time domain resource shared by the UEs without blocking 
them. 

For each BS in the under-loaded tier k, we assume that the 
location of the used bandwidth bLik in band i is uniformly and 
independently selected from band i. Equivalently, frequency 
hopping can be assumed. We further define 



B, 



1 ,Vfc e £,Vi e M, 



(2) 



and refer to 6i t k as the load of tier k in band i. 

This bandwidth sharing model helps to capture the fact that 
the load in HetNets is often unbalanced. In particular, small 
cells have limited coverage and thus the mean number of 
UEs supported by per small cell is small, yielding light load 
&i,k << 1. In contrast, the coverage of macrocells is much 
larger and thus the mean number L^k of UEs supported by 
per macrocell is large, yielding heavy load O^k w 1. 

E. Performance metrics 

It is clear that the performance of the HetNet depends on 
how the available bands are deployed. So we introduce binary 
decision variables defined as 



1 if band i is used by tier k; 
otherwise, 



for alH G M , k E K,. We group safe's into a vector x denoting 
the band deployment configuration. 

Given the band deployment configuration x, the signal to 
interference-plus-noise ratio (SINR) in band i provided by tier 
k can be computed as 



SINR,,, 



i,ko 



Yk, 



-Wi 



where Wi — bNo/Ci with Nq being the power spectral density 
of the background noise, and <ti k is the PPP of density 9i k\k 
thinned from $k- 

We are now in a position to introduce the performance 
metric used in this paper. Rate is the paramount metric in CA 
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which aims to provide UEs with very high data rate. So the 
metric used is the UE ergodic rate which measures the long 
term data rate attained by a typical UE. We next explicitly 
define the UE ergodic rate for multi-flow CA and single flow 
CA, respectively. 

1) Multi-flow CA: Using adaptive transmission schemes 
and treating interference as noise, the multi-flow UE ergodic 
rate (nats/sec) equals 



(3) 



R = X) Yl 4 K,fc^[log(l + SINR i!fc ) 
Here pi k denotes the multi-flow admission probability of tier k 

(i) 

in band i. -k, ' is the probability of the event that the typical UE 
connects to tier k in band i. The following conditions need to 
be satisfied for proper association distributions: ~ 



l,Vi, and nu > 0,Vfc, i. Note that 7r£ also represents the 
fraction of UEs served by tier k in band i and thus we refer to 
it as coveragJ^. In summary, the UE ergodic rate is determined 

(i) 

by three terms: admission probability pi tk , coverage ttY and 
spectral efficiency E[log(l + SINR^.)]. 

2) Single flow CA: As in multi-flow CA, the UE ergodic 
rate of single flow can be computed as 



R =zZ %k Yl P 4 ,fcME[log(l + SINR iifc )] 



(4) 



Here pi jk denotes the single flow admission probability of 
tier k in band i. ir k is the probability of the event that the 
typical UE connects to tier k. As a proper probability measure, 
association distribution {iTk}keK, satisfies that 2~2keSC Wk = ^' 
and ir k > 0, Vfc. 

The difference in UE ergodic rate for multi-flow and single 
flow results from the difference in the association policies. 
Recall that multi-flow UE can connect to multiple tiers simul- 
taneously and performs CA correspondingly. So the typical 
UE can pick the "best" tier in each band to connect to. 
In contrast, single flow UE can only connect to one single 
tier and performs CA with respect to the bands used by the 
selected tier. Though the former is more advanced, it may 
not necessarily lead to better performance, which is a bit 
counter-intuitive. The main reason is that better received signal 
power is not equivalent to higher UE ergodic rate which also 
depends on many other factors. We will explain this fact more 
in Section |Vl 

III. Multi-flow analysis 

In this section, we first derive formula for the multi-flow 
association distributions and admission probabilities. We then 
derive analytical expression for multi-flow UE ergodic rate. 
Finally, we study the corresponding lower bound and perform 
mean value analysis. 

A. Association distributions and admission probabilities 

Before deriving the association distributions, we first make 
the following assumption. 

2 The coverage defined here should not be confused with the usual definition 
of coverage probability, which equals one minus outage probability. 



Assumption 1. Each band is used by at least one tier and 



E 



keK, X i,k 



> 1, Vi, and 



each tier uses at least one band, i.e., 

y... m > i,vfe. 

There is no loss of generality in the above assumption since 
one can simply exclude the unused band and/or the useless tier. 
Under this assumption, the multi-flow steady state association 
distributions characterizing the network coverage are given in 
Lemma Q] 

Lemma 1 (Multi-flow coverage). For multi-How CA, let 



M) - 



= i n k )keK.>Vi £ M., where n^' is the fraction of UEs 



(i) 



that connect to tier k in band i. Then, 



g^-y • Xi,k^k{ZkPi,k) 



,Vfc e K,Vi G M, (5) 



is the normalization 



where G w = J2ieK x i,lXl( z l p i,l . 
constant. 

Proof: See Appendix iBl ■ 
The association distributions are spatial averages and are 
also the tier connection probabilities of the typical UE. In other 
words, 7Tfc also denotes the probability that the typical UE is 
"covered" by tier k in band i. Since on > 2, the association 
distributions are more sensitive to the variation in the BS 
densities than to the variation in the transmitter powers, which 
is also true for the throughput in wireless packet networks |20|. 
We further remark that Lemma[T]can also be derived following 
the approach used in lfl4l . However, the proof here may be 
more intuitive. 

With Lemma Q] we next derive the admission probability 
of each tier in each band. To this end, we assume that if a 
BS of tier k is fully-loaded in band i, then it uniformly and 
independently picks those UEs in its coverage to block. 

Lemma 2 (Multi-flow admission probability and load). In 

multi-flow CA, the admission probability and load (pi fc, #i ft) 
of tier k in band i are given by 



(Pi 



1- 



if < Di. k < & 



,6A<")(Z fc P 1 . fc ) T 

(0,0) 



-,1 ifD, k >^; 



otherwise, 



where D a = ■ x^ k \^{Z k P ltk ) £ . 

Proof: Let LT^ (x) denote the probability that the UE at 
position x connects tier k in band i. Then the mean number 
of UEs attempting to connect to a BS of tier k is given by 



Jr 2 ni (x)$w (dx) _ r o °° n« (r, e)\( i 



rdrdO 



p2tt 



Jo" 7(T ' Xkrdrde 



X k ^ Jo" rdrdO G« 



•Xi, k \M(Z k P itk y> 



where the third equality follows from the ergodicity of PPP 



,(0 



(i) 

Vx, and plugging ir k yields the last 



< 4 1 , tier k is 



i.e., n^ j (x) 
equality. So if < ■ x i)k X^ u \Z k Pi 
under-loaded in band i and thus the corresponding admission 
probability equals 1. If -^jj ■ Xi tk \^ u \Z k Pi tk )~ > tier k 
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is fully-loaded in band i and the associated UEs are randomly 
and independently blocked and thus the corresponding admis- 
sion probability equals — Similarly, we can 

derive the load 0^. ■ 
Lemma [2] describes the intuitive impact of biasing on 
admission probability and load. If a tier is lightly loaded, 
increasing its biasing factor helps to increase its load and 
thus makes the tier accomplish more. If a tier is fully-loaded, 
decreasing its biasing factor helps to increase the admission 
probability of those UEs in its coverage and thus improves 
UE throughput in the long run. 

B. UE ergodic rate 

In this subsection, we derive the analytical formula for 
multi-flow UE ergodic rate defined in (01. 

Proposition 1 (Multi-flow UE ergodic rate). The multi-How 
UE ergodic rate is given by 

/■OO rOO 

R=^2J2 2nb Pi^k / g^\t,r)dtrdr (6) 



where 



_ I .°i( e t_l) P r fc 1 « ;i -r 2 -7r(Z fc P i , fc ) «i (G« +G[ i} (t)) 



with Gfit) = E Ie /c0i,jA, {ZiPi,i)^ p{Z k Z^ 1 {e t - l),a ( ) 
andp(y,ai) — y~ J^ji. — Kq-ds. 

V <*« l + S 2 

Prop. Q] explicitly characterize the impact of band deploy- 
ment configuration x and biasing factors Z on the multi-flow 
UE ergodic rate. To obtain some insights, we next focus on the 
special interference-limited scenario where noise is ignored, 
i.e., No = 0. Under this assumption, the following Corollary 
[T]can be derived from Prop. Q] 

Corollary 1 (Interference-limited multi-flow UE ergodic 
rate). If No = 0, multi-How UE ergodic rate is given by 

roo fit) 

fl=EEft.#/ • • ■ '"■ (7) 



Corollary Q] gives a clean characterization of the impact 
of biasing on UE ergodic rate. Consider band i and assume 
there are only tier k and I, representing small and macro 
cells, respectively. The UE ergodic rate is determined by three 
terms: admission probability p^k, coverage tt^. ' and spectral 



efficiency 



G< 



GW+GW(t) 



dt. Increasing ZkZ, 1 moves more 



(i) 

UEs into the coverage of small cells, i.e., n^. ; increases. 
Meanwhile, admission probability p.^k still equals 1 since the 
small cells are typically under-loaded. The spectral efficiency 
of small cells decreases with increased Z^ZT , The reason for 
the loss of spectral efficiency is that BSs in small cells need to 
serve UEs which are moved from macrocells by biasing and 
thus generally have longer radio links than the original UEs 
in the small cells. 

In contrast, as Z^Z l increases, coverage ' of the macro- 
cells decreases, while admission probability pn of macrocells, 
which is less than 1 since the macrocells are often fully-loaded, 



increases because fewer UEs attempt to connect to macrocells. 
However, with increased ZkZ^ 1 , the change of the spectral 
efficiency of macrocells depends on specific network param- 
eters and is uncertain. This uncertainty stems from the fact 
that small BSs are uniformly distributed within the macrocell. 
Thus, UEs moved into small cells by biasing can have either 
long or short macro radio links, and correspondingly variable 
interference levels. 



C. Lower bound and mean value analysis 

In this subsection, we derive a lower bound and then based 
on mean value analysis give an approximate but much simpler 
expression for UE ergodic rate. To obtain the desired lower 
bound, we first note that the integration limits of p(y,ai) are 

2_ 

from y to co since the closest interferer is at least at 
some distance in cellular networks conditional on connection. 
If we ignore this and instead let the integration limits be 
from to oo, we overestimate the interference and thus get a 
conservative result. In particular, 

1 



p(y,on) < y 7 



ds = 2y 

1 + ST ' Jo 1 + U a * 



du 



2wy~ 



— = p(y,oti). 



(8) 



Using this inequality, the following corollary follows immedi- 
ately. 

Corollary 2 (Lower bound on multi-flow UE ergodic rate). 

A lower bound Rib on multi-How UE ergodic rate can be 
obtained by replacing p{y, a, t ) used in gf^} (t, r) with p(y, ai). 

Next we adopt the following mean value approximation to 
derive a simpler expression for UE ergodic rate: 



i£M keK. 

i^i^wm .ill*-*]- ,, m 

E[Ii\Ki = k\+Wi 

where Ki is the random tier that the typical UE connects to 
in band i, and 

h = J2 x ^ J2 p i,i H i,Y\\ Y \\~ ai ■ 

leK Keii.iYn.0 

From the proof in Appendix ICl we know the distance || Yk \\ 
conditional on Ki = k is Rayleigh distributed as 

P(||l*oll < = k ) = l-e- nr ^ ZkPi ^~^^i^ x ^ z ' p ^ 
with the moments 

r(i + ?) 



Substituting m = 1 in the above equation yields the condi- 
tional mean distance 



EiWYkJlK, = k] 




(10) 
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' ^- is the mean distance to the closest BS of tier 



Note that |, 

k from the typical UE without conditioning. Interestingly, the 
typical radio link in tier k conditional on Kj — k equals this 



distance |- 



multiplied by a factor of \l tvY 1 that is less 



than 1. In other words, conditioning reduces the mean distance 
to the closest BS, which is intuitive. 

Unfortunately, an exact computation for the conditional 
mean interference E[Ij|i<Q = k] is tedious. Instead, we 
introduce a fitting parameter 8 and approximate E[/i|i£"j = k] 
as follows. 



E[h\Ki=k] 



Y,^x i>l 6 i>l \ l P itl -E[H iiY ]- I 



~ a *dr 



**K[||Y- feo |||.K- i =ft] 



2tt 



4A fc 



(ID 



where we can see the fitting parameter 8 is used to tune the 
guard region. In particular, only those BSs outside the guard 
region of radius <5*_E[||Yfc [[ |JQ = fc] contribute to the cochan- 
nel conditional mean interference. Finally, plugging ( TTOb and 
(fTTT i into (0 yields the following approximate expression for 
multi-flow utilization factor: 



*= E E4 4 W-Mo g {i + z,7 i (2yGwr- 



(- 



2tt 



4A* 



' ^XirfijXiPij + Wi) 1 }. 
leic 



(12) 



Compared to the counterpart R given in Prop. Q] the 
approximate expression R derived from mean value analysis 
has a much simpler form which involves no integration and 
may be more friendly to use, e.g., it may serve as the objective 
function for various network optimization purposes. However, 
the fitting parameter 8 is a prior unknown. Thus, when 
applying the approximate expression R, appropriate fitting 
parameter 8 needs to be determined for example by comparing 
the approximate expression to the exact one. A detailed study 
on selecting the fitting parameter 8 is beyond the scope of this 
paper. 



IV. Single flow analysis 

In this section, we derive formulas for the single flow 
association distribution and admission probability and UE 
ergodic rate. The analysis is similar to that of multi-flow and 
thus this section is brief. 

For single flow, we make the following additional assump- 
tion for tractability. 

Assumption 2. In single flow CA, the transmit power of each 
tier k is not band-dependent, i.e., P^k = Pk, Vi 6 M. 

This assumption is slightly restrictive and may be relaxed 
in future research. The single flow association distribution is 
then given in Lemma [3] 

Lemma 3 (Single flow coverage). For single flow CA, let-rr = 



{^k)keK.M where 7r fc is the fraction of UEs that connect to tier 
k. Then, 



7r fc = 27rAfeGfc, Vfc e K., 



(13) 



where Gk — / °° rhk(r)dr with 



Z,Pi . _a_ 



h k (r) = exp(~l^J2 x l(lT L T^) a, '■ r ° ! * ) 

ie/c k k 

andk* = aigmm l£ M:x, :k ^o cti . 

Proof: See Appendix ID] ■ 

With Lemma [3] we next derive the admission probability 
of each tier which uniformly and independently blocks UEs 
in its coverage if it is fully-loaded. The proof of Lemma [4] is 
similar to that of Lemma [2] and is omitted. 

Lemma 4 (Single flow admission probability and load). In 

single flow CA, the admission probability and load (pi t ki 9%,k) 
of tier k in band i are given by 



11 ,1) if2ir\^Gk>^ 



(0, 0) ifx hk = 



Finally, we derive the analytical formula for single flow UE 
ergodic rate defined in (0|. 

Proposition 2 (Single flow UE ergodic rate). The single flow 
UE ergodic rate is given by 

fi=Vr2* iPi J / g i l s l(t,r)h k (r)rdrdt, 



ieM keK. 



(14) 



where 



; } (t,r) = exp{-r a *(e* - l)P fe "V - 7r^0 M A, 



leK, 



ZtP 2 (Z k . t 



on) -r^ }. (15) 



K z k Pk' ' {zr 



Proof: As suggested in Appendix lAl we need to find the 
radius r k .i of the guard region and the single flow distance 
distribution f\\Y ko \\( r \ K = k) of the associated BS Y ko in 
tier k. The former can be determined from that Z k Pk-r~ ak * = 
ZiPi ■ r^" L * . As for the latter, note 



P(||n || >r\K = k) 



= — 27rAfc / ./■ • ( 

TTfc 



P(||n [|>r,A' = fc) 



dx, 



where the last equality can be shown by using similar argu- 
ments as in the proof of Lemma [4] Hence, f\\Y ko \\( r I K = 
k) = -^-2nXkrhk(r),r > 0. The remaining proof follows as 
in the multi-flow case. ■ 



3 We use 7r to denote both the single flow association distribution and the 
constant Pi. The meaning should be clear from the context. 
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TABLE I: Simulation parameters 



Fixed UE density with X iu) =24 



Density of macrocells Ai 


(7T0UU ) 


Ti£»ncitw c m oil j^ollc \ 
UCllalLV Ul allld.ll OClla A 2 




n^ncit\7 r»f TTPc \{ u ) 


9/i v fVcpin 2 "!— i 


iviax 1a power 01 macro dos 


w 


Max Tx power of small BSs 


1 w 


INOlSe rol_J 


1 7/MRm 
— 1 / 4Ur>m 


800MHz band's wavelength [i\ 


0.375m 


800MHz band's path loss exponent ai 


3 


OUUIVIOZ UallLl a UallLlW 1LIL11 ±J\ 


QIVTH7 


2.5GHz band's wavelength /i 2 


0.125m 


2.5GHz band's path loss exponent a 2 


4 


2.5GHz band's bandwidth B 2 


9MHz 


UE min bandwidth requirement b 


1.8MHz 


Band deployment 


[1,1; 1,1] 



V. Numerical results 

In this section, we provide some numerical results to demon- 
strate the analytical results. The set-up is a typical HetNet 
consisting of 2 tiers (macro and small cells) and 2 bands 
(800MHz and 2.5GHz band). The specific parameters used are 
summarized in Table U unless otherwise specified. We further 
introduce a shorthand notation [xi.i, £2,1; £1,2 , £2, 2] to denote 
the band deployment configuration for ease of description. For 
example, configuration [1, 1; 1, 0] denotes that tier 1 uses both 
band 1 and band 2 while tier 2 only uses band 1. 



A. The need for biasing 

In this subsection, we show that, though small cells can 
bring significant capacity increase to the cellular networks, bi- 
asing can make small cells accomplish more if the UE density 
is not large. To explore this issue, we show the network sum 
rate, which equals UE ergodic rate multiplied by UE density, 
as a function of normalized density of small cells in Fig. [3] In 
the top subfigure we fix UE density as increasing the density 
of small cells, while in the bottom subfigure UE density 
increases proportionally with the density of small cells. Both 
two subfigures demonstrate that network capacity increases 
almost linearly as the density of small cells increases, which 
shows the promise of small cells. Furthermore, the following 
observations can be made. 

1) : In the top subfigure small cells become more and more 
lightly loaded because we blindly increase the density of small 
cells but the UE density does not grow correspondingly. In 
this case, we can see that biasing helps increase sum rate. For 
example, when A2/A1 — 4, the sum rate attained with 30dB 
biasing towards small cells is about 35% larger than the sum 
rate attained without biasing. 

2) : In the bottom subfigure the proportional increase of UE 
density makes small cells already have heavy user load when 
A2/A1 > 4. In this case, we can see that biasing does not help. 
Worse still, biasing reduces the sum rate when the small cells 
already have heavy user load. This implies that using a load 
balanced model to study biasing is inappropriate. 




T\lormalized densityof small cells: 



Varied UE density with X (u) /X z = 12 




ormalized density of small cells: X JX 



Fig. 3: Impact of the density of small cells on network sum 
rate: Biasing factor Z\ of macrocells is kept constant, i.e., 
Z x = OdB. 



3): Comparing the two subfigures we see that sum rates 
attained in the top subfigure is significantly smaller than their 
counterparts in bottom subfigure. This is because small cells 
are under-utilized in the top subfigure, which results in low 
sum rate. This observation gives us the intuitive design rule 
that adding small cells should be demand-based and too many 
small cells will lead to waste. 

B. The complex impact of biasing 

In the previous subsection, we showed that biasing can 
help small cells accomplish more if the network load is 
unbalanced. In this subsection we provide more numerical 
results to carefully study the impact of biasing on the HetNet. 
For brevity, we first assume only 800MHz band exists and 
is used by both macro and small cells. In this singe-band 
setting, multi-flow CA is the same as single flow CA. The 
impact of biasing on coverage probability, load and admission 
probability is shown in Fig. H] It can seen that the typical 
impact of biasing on these parameters is consistent with our 
previous high level discussions in Section IIII-BI when Z 2 is 
"typical" in the sense that Z 2 £ [0dB,21dB]. Note that, as 
shown in Fig. [4] the load of small cells is 1 when Z 2 = lOdB. 
However, this by no means implies that lOdB biasing gain is 
enough. In particular, the network still can still benefit from 
further increasing biasing factor Z 2 to shift more UEs from 
macrocells to small cells. 



8 



The impact of biasing on spectral efficiency is shown 
in Fig. [5] where the situation becomes more complicated. 
Consider small cells first. When Z2 is typical in the sense that 
Zi e [0dB,21dB], the top subfigure shows that the spectral 
efficiency of small cells decreases as Z2 increases. This is 
because small cells need to serve UEs of longer radio links 
due to biasing, as explained in Section IIII-BI Interestingly, 
increasing Z2 further beyond 21dB leads to increased spectral 
efficiency of small cells. This is because macrocells become 
more lightly loaded as Z2 increases beyond 21dB (ref. Fig. 
3J. As a result, small cells experience less interference from 
macrocells, which outweighs the loss resulted from serving 
UEs of longer radio links. In contrast, the bottom subfigure 
shows that the spectral efficiency of small cells decreases as 
Zi increases. We do not observe that the spectral efficiency of 
small cells increases after Z2 increases beyond some threshold 
value. This implies that for interference-limited small cells, 
as Z 2 increases, the gain from reduced cross-tier interference 
does not outweigh the loss from serving UEs of longer radio 
links. 

The conserve is true when it comes to the spectral effi- 
ciency of macrocells, as shown in Fig. [5] In particular, when 
Z-2 < 21dB, the top subfigure of Fig. shows that the 
spectral efficiency of macrocells increases as Z2 increases, 
while the bottom subfigure shows that the spectral efficiency of 
macrocells first decreases and then increases as Z2 increases. 
The observations here are consistent with our discussion in 
Section IIII-BI Interestingly, as Z2 increases beyond some 
threshold value (i.e., lOdB in Fig.[5]l, the spectral efficiency of 
macrocells keeps increasing as Z2 increases. This is because 
only a few UEs of good geometry are served by macrocells. 
Nevertheless, this comes at the cost of the spectral efficiency 
of small cells. 

Next we show the overall impact of biasing on UE ergodic 
rate. For brevity we focus on multi-flow CA in Fig. [6] As the 
biasing factor Z2 of small cells increases, the UE ergodic rate 
first increases and then keeps almost constant for a certain 
range of Z2 and finally decreases. This can be understood 
by examining the detailed explanations about the impact 
of biasing on the various aspects of the HetNet including 
coverage probability, load, admission probability and spectral 
efficiency, which together determine the UE ergodic rate. It is 
observed that the maximum biasing gain is about 25% to 40%, 
attained at large biasing factors (ranging from 15dB to 30dB). 
This demonstrates that we need to work on large biasing factor 
for small cells in real HetNets to reap the benefit of range 
expansion. 



C. Band deployment 

In this subsection, we provide numerical results to study 
how different band deployment configurations affect the Het- 
Net performance. For brevity, we first focus on single flow CA 
and show the impact of different band deployment configura- 
tions on the single flow coverage 7r, load 9 and admission 



800MHz band 




10 15 20 25 30 35 40 






— b— Macrocells 




— *— Small cells 




Biasing of small cells Z 2 (dB) 

Fig. 4: Impact of biasing on coverage probability, load and 
admission probability. 

800MHz band 




Biasing of small cells Z 2 (dB) 
800MHz band: Interference-limited 




Biasing of small cells Z 2 (dB) 



Fig. 5: Impact of biasing on spectral efficiency. 



probability p in Table iHn We remark on the numerical results 
in Table [TT] as follows. 

Deploy the 800MHz band in macrocells to ensure coverage. 
Examining the coverage under band deployment [0, 1; 1, 0], we 
find that the coverage probability of the macrocells is only 
0.14. This implies that in our model 86% of the UEs have 
to connect to the small cells. While feasible in theory, in 

4 Note that by definition 9i k = 82 ki k = 1> 2 unless either X\ k or x 2 k 
is zero. For simplicity, in Table [TT] we only show the nonzero value which 
represents both 6\k an d ^2 fe> while by default we understand that 0f k is 
if x i,k = 0. Similar rules apply to Pi,k- 



9 





Deployment 


[0,1; 1,0] 


[1,0; 0,1] 


[1,0; 1,1] 


[1,1; 1,0] 


[1,1; 0,1] 


[1,1; 1,1] 


Coverage 


Macrocells 


0.14 


0.99 


0.85 


0.85 


0.99 


0.85 


Small cells 


0.86 


0.01 


0.15 


0.15 


0.01 


0.15 


Load 


Macrocells 


0.66 


1.00 


1.00 


1.00 


1.00 


1.00 


Small cells 


1.00 


0.03 


0.35 


0.35 


0.03 


0.35 


Admission 


Macrocells 


1.00 


0.21 


0.24 


0.24 


0.21 


0.24 


Small cells 


0.48 


1.00 


1.00 


1.00 


1.00 


1.00 



TABLE II: Impact of band deployment on single flow coverage, load and admission probability 




5 10 15 20 25 30 35 40 



Biasing of small cells Z 2 (dB) 

Fig. 6: Impact of biasing on multi-flow UE ergodic rate and 
utilization factor. 



reality the 800MHz band should be deployed in macrocells 
to ensure coverage. Thus, in the sequel we only study those 
band deployment configurations with x\\ = 1. 

Deploy the 800MHz band in small cells to better balance 
load. Comparing the deployment [1,0; 0,1] to [1,0; 1,1], we 
can see that the coverage and load of small cells are very 
low if the 800MHz band is not used in the small cells. 
So we suggest that the 800MHz band be deployed in small 
cells too. Similar conclusions can be drawn by comparing the 
deployment [1,1; 1,0] to [1,1; 0,1]. 

We next show the impact of band deployment on single 
flow spectral efficiency in Table iHll Generally speaking, band 
deployment affects the spectral efficiency through its influence 
on the distribution of the network load, which further affects 
the desired signal power as well as interference. For example, 
let us compare the deployment [1,0; 0,1] to [1,0; 1,1]. With 
deployment [1,0; 0,1], the load of the small cells is very 
low (as shown in Table |ll|i and thus they only need to serve 
those UEs which are very close to them. Besides, since the 
macrocells do not use the 2.5GHz band, there is no cross-tier 
interference. These two factors lead to the very good average 
spectral efficiency of the small cells in the 2.5GHz band, i.e., 
8.23 bps/Hz. In contrast, with deployment [1,0; 1,1], small 
cells have much higher load and need to serve UEs with longer 
radio links, resulting in lower spectral efficiency: 3.26 bps/Hz. 
Comparing further band deployment [1, 0; 1, 1] to [1, 1; 1, 1], 
we observe an even lower spectral efficiency (i.e., 2.4 bps/Hz) 



of small cells in the 2.5GHz band due to the additional cross- 
tier interference from the macrocells. 

Now let us study the overall impact of band deployment on 
single flow UE ergodic rate in Table llVl which lists all the R j 
as well as their sum R. We observe that the UE ergodic rate 
is highest with full deployment [1, 1; 1, 1], while it is lowest 
with orthogonal deployment [1, 0; 0, 1]. We further notice that 
orthogonal deployment [1, 0; 0, 1] has lower sum rate than the 
rate 1.57Mbps (not shown in Table II Vb attained by deployment 
[1,1;0,0]H which corresponds to a 1-tier cellular network 
equipped with the CA feature. This demonstrates that splitting 
the spectrum between different tiers in HetNets is not a good 
strategy. This insight is consistent with that observed by 
industry. Indeed, HetNets can even have worse performance 
under spectrum splitting between macrocells and small cells 
as reported by Ericsson [8|. The reason is that there is a 
capacity loss in macrocells by spectrum splitting and this loss 
can only be partially overcome by adding a large amount 
of small cells. Thus, HetNets should allow spectrum reuse 
augmented with appropriate intercell interference management 
ifTUl . Also, comparing the rates provided by small cells under 
different band deployment configurations, we see again that 
the small cells do not accomplish much because of their 
limited coverage if they do not use the 800MHz band. 

D. Multi-flow vs. Single flow 

In this subsection, we compare the performance of multi- 
flow to that of single flow under different band deployment 
configurations. Fig.|7]shows the multi-flow and single flow UE 
ergodic rate as a function of the density of small cells. We can 
see that [1,0; 1, 1] is almost as good as [1, 1; 1, 1]. However, 
[1,0; 1, 1] has bad performance in multi-flow CA in the low 
small cell density regime, i.e., A2/A1 < 5, while [1, 1; 1, 1] 
is optimal or near optimal in both multi-flow CA and single 
flow CA. This suggests that both the macro and small cells 
use both the 800MHz and 2.5GHz band. However, most of 
the current macro BSs only support transmissions in 800MHz 
band. If so, we further suggest that [1,0; 1,1] be used. This 
would still provides very good performance because most of 
the current UEs are single flow based and thus the bad multi- 
flow performance of [1,0; 1, 1] in the low small cell density 
regime would not bother much, at least for the time being. 

Comparing the two subfigures in Fig. [7] we observe the 
interesting fact that multi-flow CA does not necessarily provide 

5 Though [1, 1;0, 0] does not satisfy our Assumption [TJ in Section ITTT1 our 
results can be readily applied to this scenario by simply ignoring the small 
cells. 



10 





Deployment 


[1,0:0,11 


[1,0; 1, 11 


[1, 1; 1,01 


[1,1; 0,11 


[l, l; 1, il 

L ) 5 5 J 


Macrocells 


800MHz 


1.30 


1.35 


1.35 


1.30 


1.35 


2.5GHz 


0.00 


0.00 


1.75 


1.62 


1.69 


Small cells 


800MHz 


0.00 


1.35 


1.35 


0.00 


1.35 


2.5GHz 


8.23 


3.26 


0.00 


6.91 


2.50 


TABLE III: Impact of band deployment on single flow spectral efficiency (bps/Hz) 




Deployment 


[1,0; 0,1] 


[1,0; 1,1] 


[1,1; 1,0] 


[1,1; 0,1] 


[i,i; 1,1] 


Macrocells 


800MHz 


0.70 


0.73 


0.73 


0.70 


0.73 


2.5GHz 


0.00 


0.00 


0.95 


0.87 


0.91 


Small cells 


800MHz 


0.00 


0.51 


0.51 


0.00 


0.51 


2.5GHz 


0.27 


1.23 


0.00 


0.23 


0.91 


Sum 




0.98 


2.58 


2.19 


1.80 


3.06 



TABLE IV: Impact of band deployment on single flow UE ergodic rate (Mbps) 



Multi-flow CA 




5 10 15 20 25 

Normalized density of small cells: ^JK 
Single-flow CA 




Fig. 7: 



5 10 15 20 25 30 

Normalized density of small cells: ^JK 

Impact of band deployment on single flow UE 
ergodic rate 



better performance than single flow CA. This might seem a 
bit counter-intuitive. We use the deployment [1,1; 1,0] as an 
example to explain this observation. With multi-flow CA, all 
the UEs in the network attempt to connect to the macro tier 
in the 2.5GHz band, resulting in only a small share of the 
resource for each UE. Besides, many of these UEs have long 
radio links and correspondingly weak signal power. So we can 
improve the UE ergodic rate if we could prevent part of these 
UEs from connecting to the macro tier, which is implicitly 
implemented in the single flow CA. In particular, those single 
flow UEs connecting to the small cells would not pursue to use 



the 2.5GHz band. Though these UEs suffer from a bit rate loss 
without accessing the 2.5GHz band, this significantly benefits 
other UEs connecting to the macrocells since they generally 
are UEs of good geometry w.r.t. to the macro tier and now can 
get a much larger share of the resource in the 2.5GHz band 
and thus enjoy much higher rate. We conclude that this overall 
results in better performance for single flow CA in terms of 
UE ergodic rate. 

Before ending this subsection, we remark that Fig. [7] focuses 
on a specific scenario: no biasing and interference mitiga- 
tion. In reality, deployment choices need not be restricted to 
[1, 1; 1, 1] or [1, 0; 1, 1] and there will be many operating points 
in between. 

VI. Conclusions 

This paper focuses on how to best exploit carrier aggrega- 
tion in HetNets. To this end, we have proposed a general yet 
tractable M-band K-tier model for carrier aggregation which 
captures the key network parameters. In particular, this model 
incorporates an appropriate notion of load which allows the 
study of biasing, as well as the impact of band deployment. 
Under this model, we derive analytical expressions which 
explicitly express the dependency of the network metrics 
on the biasing factors and band deployment configuration. 
This work can be extended in a number of ways. In par- 
ticular, many other distinct features such as multiple-input 
multiple-output (MIMO) technique, coordinated multi-point 
transmission/reception (CoMP) and device-to-device (D2D) 
communications should be jointly studied with CA. 

Acknowledgment 

The authors thank Rapeepat Ratasuk and Bishwarup Mondal 
of Nokia Siemens Networks for their help in identifying the 
problem and refining system model and numerical results. 

Appendix 

A. The CCDF of SINR of band i in tier k 

We derive the CCDF of SINR of band i in tier k, which lays 
the foundation for the analysis in this paper. The derivation is 
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mainly based on the tool of stochastic geometry and follows 
the main technique used in fl4l . [18|. 



To this end, we first write the CCDF of SINR of band i in 
tier k as 



P(SINR a > T) = 

,Pi,kH itko r~ 



E 



'(- 



>r) 



fm, k0 \\(r)dr, 



(16) 



1% + Wi 

where h = Y^ieK. 1 *! with hi = T,Ye& it i\Y ho p U H i,y I 
Y ||~ Qi , and f\\Y ikQ \\{ r ) is the distance distribution of the 
associated BS Yi :ko . Note kQ \\(r) depends on the specific 
association strategies. We will derive f\\Yi ko \\( r ) later for 
multi-flow CA and single flow CA, respectively. Here, 



E 



Pi.kHi ka r~ 



>r) 



= E [P(H lM > r^TPr^Ii + Wi )) 

= E[cxp{-r a 'TPr k 1 (I l +w l )) 

= eM-r a ^P^ Wl ) J] £ it i(r ai TPr£) 



(17) 
(18) 
(19) 



where A,z(s) = E[e is the Laplace transform of 

h,l> ( TT81 follows from the fact that Hi yko ~ Exp(l), and 
(fT9l follows from the independence of the K underlying 
PPPs. Following similar arguments as in [18], we have that 

d,i(r a *TPr^) equals 



E 



'^■^«„A % Pi,iH tlY \\Y\\- 



n 



E \ e - rC " r ^-P^\\Y\\~ a 'H t , 



=E 



n 

aire*,, a, / r °° 



1 

i + r«-rp 4 ,p. 



-<// 



(20) 



(21) 



(22) 



(23) 



In (fJTJ and (l22l i. we use the independent and identically 
Rayleigh distributed assumption on the interference fading, 
i.e., H i>Y ~ Exp(l), and the fact E[e~ sH ] = ^ where 
H »j Exp(A). In ( f23l , we apply the probability generating 
functional of PPP (see, e.g., ifTBI ). Note the last integration 
limits are from r^h,l to oo since, conditional on association 
with some tier k in band i, the closest interferer in tier I is 
at least at some distance r^h,i which depends on the specific 
association strategy. We will specify r^i later. Now ( fT9b can 
be further simplified as 



E 



P( 



Pi,kHi,k i 



-r a iTP7 



>r) 



(24) 



Furthermore, 



t 



< 



. ^ ! + r - air -xp-i p . ik . tl 
t 



-dt 



o 1 + r-^T-iprip^ ■ t<** 
Thus, we can lower bound (124-b as 

Pi, k H iM r- a 



at = -i— — • r . 



ohm*) 



E 



>r) 



> e 



2 «- «»,i^i( r n,i p < fc > 



(25) 



If we further obtain fco || (r) and r^fc^, we can evaluate the 
CCDF of SINR of band i in tier k (resp. its lower bound) by 
plugging ([24j (resp. (|25])) into ( [TBI 

5. Proof of Lemma Q] 

The association distribution is equivalent to the tier connec- 
tion distribution of the typical UE. Let Ki be the random tier 
that the typical UE connects to in band i. We next compute 
its probability mass function (pmf) as follows: 



.(0 



\Ki = k) 
= P(Z k P^ k \\Y ko \\ 



max^P^Hrj-^)) 
leic ' 



= P(ir\\Y ko \\ 2 (Z k P hk ) ^ =mln(ir\\Y l0 \\ 2 (Z l P itl )-^)) 

It is known that in PPP || Yj 1| is Rayleigh distributed [21 ], i.e., 
the probability density function (pdf) of \\Yi \\ is given by 

f llYlJ (r)=2n\ l re- x ' r2 , r > 0, 

from which we can derive that 7r||Yj || 2 (^Pij) ~ ~ 
Exp((ZiPij)^i Xi), yi G /C. These exponentially distributed 
random variables are independent since the underlying K PPPs 
are assumed to be independent. It follows that the probability 
that the fc-th one is the smallest among the K independent 
exponential random variables is 



\K % = k) 



\k{Z k Pi, 



(26) 



Taking into account the configuration of the band deployment, 
we have 



4° = F(Kt = k) = -L . x i!k \ k {Z k P itk )% , 



(27) 



where the denominator cannot equal due to our assumption 
that each band is used by at least one tier. 

C. Proof of Proposition Q] 

Denote Ri\ Kz =k = E[log(l + SINR^)^ = k]. Note that 
E[log(l + SINRi, fc )|tfi =k] 

I'.jU.j r 1 



log(l 



fWY^wiAKi = k)dr. 

(28) 
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where f\\y k \\{ r \Ki — k) can be derived as follows: 

P(\\Y ko \\>r\Ki = k) 

= F(7r\\Y ko \\ 2 (Z k P l , k y^ > 7rr 2 (Z k P itk )-^\Ki = k 



(29) 



P (wmfrllY^ffrPi,,)--* :le)C)> nr 2 {Z k p, k ) 

(30) 

[]p( 7 r||y io || 2 (Z i P iJ )-^ > nr 2 (Z k P itk )-~ ' 



(31) 



(32) 



where (|3T1 i follows since the exponential random variables in 
question are independent. Taking into account the configura- 
tion of the band deployment, the pdf of \\Y k() || conditioned on 
Ki = k is given as follows. 

f\\Y k j(r\K t = k) 

= 2irr(Z k P ltk f^G^ ■ e-* r2 ( z * p '-^~ G ' M ,r > 0. (33) 
Also, 

I'.j ll. i r " - 



E 



log(l + 



E 



h + Wi 
log(l + 



Ii + Wi 



)>t 



dt 



E 



(H itko > r a *(e t -l)Pr k 1 (I i +w i )) 



dt. (34) 



By comparing (l34l to (jT7j, we note the only difference is that 
the term T in dTTi l is replaced by e* — 1 in d34l i. So we can 
directly utilize the result given in (l24l . We further simplify 
(ED by using the fact that r i)fc ,j = (Z t Pu)^ (Z k P hk )~^ r 
in multi-flow CA and changing the variable t to s = 
(r^TPijPr^-kt 2 to get 



E 



r,jii,, r " 



>r) 



-r a irF- k 1 w 4 -r 2 -7r(2 fe F i ,fc)"' s i'6W 



(35) 



where G 



(i) 



'/, HieK d iAl( z i p i,i)^ p{TZ k Z^, ai ) and 
p(y,ai) = y~ J x ^- — ^ds. Thus, plugging (E2K d33l 
and d35j (with T substituted by e* - 1) into i? = 
EieA4 SfeGK ^ft bR i\K t =k Y ields th e desired result. 



£). Proof of Lemma [3] 



In single flow CA, the typical UE scans over all the tiers 
and bands and then connects to the tier k that provides the 
strongest biased received power in some band say fc* e M. 
Then the UE performs CA with respect to this selected tier. 



Let K be the random tier that the UE connects to. Then, 

TTfe = P(K = k) 

= P(Z k P k \\Y ko \\-^ =max(Z i P i ||^.o|r Q 0) (36) 

l,i 

= F(Z k P k \\Y ko \\- a ** = max(Z,fl||y, |[- a '* : J e £)) 
= l[P(7T(Z l P l y^\\Y l0 \\ 2 > n(Z k P k )-*F\\Y k0 \\^) 



IP 

l^k 

l=ik 



n 



p\\\Y lo t> * {ZkPky 7 \\Y k0 \\^)\m 



\ TT(ZiP) "V 

c ^{-KM^W\\Y k0 \\^)\\\Y ka \ 

^k^k 



(37) 



2ir\ k r exp(—Tr\ k r 2 )dr 



(38) 



2,TAfc / rexp(-7ry^ A;( J ^ ) a i* r a i* )dr 
Jo leK Z k P k 



where we use %{Z x P itl ) ^\\Y lo \\ 2 ~ Exp((Z ; P M ) =i \) 
(|37| | and the fact that ||Yfc || is Rayleigh distributed in 
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